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Statistical properties of Fano resonances occurring in photoabsorption to highly excited atomic 
or molecular states are derived. The situation with one open and one closed channel is analyzed 
when the classical motion of the excited complex in the closed channel is chaotic. The closed channel 
subspace is modeled by random matrix theory. The probability distribution of the Fano parameter is 
derived both for the case of time reversal symmetry (TRS) and broken time reversal symmetry. For 
the TRS case the area distribution under a resonance profile relevant for low resolution experiments 
is discussed in detail. 
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Introduction — Resonance phenomena are ubiquitous 
in highly excited quantum systems. Typical examples 
are auto-ionizing atomic resonances [0 or photodissoci- 
ation of rovibrational molecular states [^. When the 
density of resonances is high a common and natural way 
to extract information about resonance properties is to 
resort to statistical methods ^, |^. Typically as a 
consequence of the high excitation the resonant complex 
behaves chaotically from a classical point of view. It is 
then appropriate to model the excited complex by means 
of random matrix theory ||, 0, || . 

Here I address a problem frequently encountered in 
photoabsorption processes of atomic and molecular sys- 
tems: What are the statistical properties of line shapes 
when photoabsorption into a highly excited resonant 
state takes place? Since the seminal work of Fano ||] 
it is known that when at least two pathways exist for 
an atomic or molecular system to decay after having 
absorbed a photon the line shapes of resonances dif- 
fer in general considerably from the Breit-Wigner form. 
In photodissociation of molecules or photoexcitation of 
an atom to an auto-ionizing resonance interference be- 
tween the indirect and the direct decay process results 
in a Beutler-Fano resonance profile of the photoabsorp- 
tion cross section [|], |l0[ Q (abbreviated as Fano profile 
in the following). The shape of the resonant part of the 
cross section is determined by a single parameter, the 
Fano parameter q. This contribution is aimed at such a 
statistical theory of Fano resonances. 

It is worthwhile to stress that atomic and molecular 
systems are ideally suited to test the predictions on Fano 
resonances as neither Coulomb or dephasing effects which 
may become relevant in mesoscopic structures have to be 
taken into account . 

Theory — For photodissociation of a molecule Fig. |l| 
schematically depicts the situation envisaged [|l3|. The 
molecule is coherently laser excited from the ground state 
or a low lying state |0) of energy Eq onto two electronic 
potential surfaces. The transition is accomplished by the 
dipole operator D. Electronic surface SI is the open 



channel — the energy E of the molecular complex is 
above the dissociation threshold — while the motion on 
the potential surface S2 is bound (closed channel). The 
dipole matrix element for transition to the open channel 



is given by di = (0|D|V'^''^^), where \4'e"'^') ^^le regu- 
lar continuum solution at energy E in the open channel, 
normalized in energy. When E coincides with the en- 
ergy En of a bound state \(j)n) in the channel 2 the closed 
channel carries the transition amplitude (O|Z?|0„) in the 
absence of coupling to the continuum. In the following 
it is assumed that transitions between the two excited 
manifolds are possible. In the diabatic representation 
the coupling between the two surfaces is given by a non- 
diagonal potential V. Then the eigenstate \(f>n) turns into 
a resonance with width r„ = 27r|(?/'lj''^-'|V^|0„)P- 

In the regime r„ <C A of isolated resonances (A be- 
ing the mean energy spacing of resonances) the oscillator 
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FIG. 1: Schematic sketch of the photodissociation process: 
Laser excitation takes place from a low lying electronic po- 
tential surface (dipole operator D). The classical dynamics 
on surface S2 is bound and chaotic. Dissociation can take 
place directly by dipole excitation to channel 1 or indirectly 
by excitation to channel 2 and the coupling V between the 
surfaces SI and S2. 
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strength Df{E) for the dipole transitfon is given by 0| 



Df{E) = Df 



(bg) 
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where Df'-^s) ^ 2^(i; - i;o)|dip/?i^ is the background 
oscillator strength for the transition to the open channel 
in the absence of coupling to the closed channel, fj, is the 
reduced mass of the excited complex and e„ = 2{E—En — 
A„)/r„ the reduced energy. The energy shift A„ of the 
resonance is neglected in the following. The line shape 
of the resonance is parameterized by the Fano parameter 
g„. For a non-degenerate continuum it is given by 0| 



+ q, q 



(2) 



The irregular continuum solution appears in the 

expression for qn because near resonance the wave func- 
tion in the open channel is a superposition of the reg- 
ular and the irregular solution. Both di and q depend 
only weakly on E and are assumed to be constant within 
the energy window over which a sample of resonances is 
taken. From (||) it is seen that the distribution of a set of 
Fano parameters {qn} taken over a stretch of the energy 
spectrum is determined by the statistical properties of 
the set of the eigenstates {!</>«)} of the closed channel. 

The distribution of the Fano parameter is derived un- 
der the two following conditions: Firstly the classical 
motion of the excited molecular complex on the elec- 
tronic surface SI of the closed channel is completely 
chaotic. This ensures that the statistical properties of 
generic wave functions in chaotic systems apply and that 
the closed channel subspace can be modeled by random 
matrix theory Secondly the excitation process and 
the coupling between the two electronic surfaces are as- 
sumed to be spatially well separated: the overlap of the 
initial wave packet (r|D|0) and the coupling potential 
V{r) in coordinate representation is negligible. There- 
fore X = {0\D\(j)n) / di and y = 7r(?/;^°^^ |T^|0„) can be 
taken statistically independent random variables [ p^ . 

Distribution of the Fano parameter — With these in- 
gredients the calculation of the probability distribution 
P{q) of the Fano parameter is straightforward. (The in- 
dex n of the resonance is omitted in the following.) In 
the case of time reversal symmetry the wave functions 
and viz. x and y can be chosen real. The closed channel 
subspace is modeled by x matrices taken from the 
Gaussian orthogonal ensemble (GOE). All results are un- 
derstood in the limit N oo where x and y are Gaussian 
random random variables with zero mean and variances 
and ay. The probability distribution is given by 
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and P^^ (x) is the probability distribution of x (and like- 
wise for y). The integral is most easily evaluated by using 
the Fourier integral representation of the delta function 
and the result is 
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-Pgoe('7; s) = - 
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(Jxloy . (4) 



The probability distribution of the Fano parameter in the 
GOE case thus turns out to be Lorentzian with mean 
value g and width s. 

The width s of the probability distribution is related 
to the coupling strength V between the closed channel 
and the continuum channel and the ratio of the dipole 
transition matrix elements to both channels. Assume V 
can be written in the form V — AVq where A character- 



izes the coupling strength and Vq is fixed. Then Oy 



A 



and therefore s ^ 1/A. For strong coupling to the con- 
tinuum the Lorentz distribution acquires a small width 
centered around g. The same holds if direct photoexci- 
tation dominates over the indirect process since then Ox 
becomes small. 

In the case of broken time reversal symmetry the 
Hamiltonian of the closed channel is modeled by x 
matrices from the Gaussian unitary ensemble (GUE). In 
this case x and y are complex Gaussian random vari- 
ables with independent real and imaginary parts and the 
Fano parameter is in general complex. Most conveniently 
the distribution of q is characterized by the probability 
distribution of the phase and the modulus of the 
quantity q— (q). The phases Lpx of x and ipy of y are 
uniformly distributed mod(27r) and the same holds for 
fq = y^x^ f y ■ Denoting the modulus of x by r^ and of y 
by ry the probability distribution PcuEifq) is given by 
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^'GUE(r,) = J dr^J dryPir,)P{ry)6l^rq - j , (5) 

where r^ (and likewise ry) has the probability distri- 
bution P{rx) = CT~^^r2. exp(— r^/2a-^^). [Notice that 
P{r^) ~ exp(— r^/2(T^ ) has the form of a Porter-Thomas 
distribution for GUE.] Proceeding as in the GOE case 
the integrals turn out to be more complicated as they 
involve error functions. The final result is nevertheless 
simple and reads 



P 



GUE 
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{rg > 0) , (6) 



where s,- = Or^l'^ry- The distribution is delta peaked 
at = and has a local maximum at rf^ — Sr/\/i 
with P(rf^) = 3V3/(16s,.). A gain, as discussed for the 
GOE case the width Sr of the probability distribution 
-Pgue (?■(?) is determined by the strength of the coupling 
between the two channels and the ratio of the strength 
between the dipole transition to the closed and to the 
open channel. 
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FIG. 2: Peak height distribution PGOE(z/g^) in the GOE 
case for different values of f3 — q/s, s being the width of 
the (/-parameter distribution Eq. (4). For /3 > 1.95517 the 
distribution displays a local maximum between zero and one. 



Peak height distribution for GOE — For the TRS case 
the maximum of the Fano profile is given at e — 1/q with 
peak height z = [Df - Df^^'s)y Dfibg) = ^2 relative to the 
background oscillator strength. The distribution of peak 
heights can easily be calculated from (H) and is given by 



GOE 



(^) 



2^ 



[i'GOE(V^;s)+FGOE(-V^;s)] (7) 



with z > 0. In Fig. g -Pgoe(2:) is plotted in terms of the 
scaled variable z /if and different values of the parameter 
/3 = q/s. It has a local maximum at a value < z/'f < 1 
for /3 > 1.95517. Therefore for a local maximum to occur 
the width s of the Lorentz profile (Q) must be smaller 
than approximately twice the average value q of the Fano 
parameter. For (3 oo the position of the maximum 
approaches z =7^ . 

Profile area distribution for GOE — In low resolu- 
tion experiments the relevant experimental quantity is 
the profile area distribution under the resonance rather 
than the peak height distribution. To be more specific as- 
sume that the band width SE of the exciting laser beam is 
larger than the width F of the resonance but still smaller 
than the mean spacing A between adjacent resonances 
(F < SE < A). Assuming a rectangular laser profile 
the excess oscillator strength averaged over the resonance 
profile with respect to the background strength is given 

by a 
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Df^''^\E')[F{q;e)~l]dE' (8) 
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where F{q; e) = \q + ep/ (1 + e^) is the profile function of 
the resonance [cf. (^]. Since F < 5E the range of inte- 
gration in (H) can be extended to infinity and D/^™'''' (E) 
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FIG. 3: Distribution P(G>) of the scaled resonance area = 
aS, see discussion after Eq. (10). The average value of the 
Fano parameter is fixed at q — and s is varied. 



is given by 



Df^""\E) ^ Df'^^^^ 
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de, (9) 



which performing the integration gives Df / Df^^^^ = 
TiT{q^ — l)/2 = 8. In the following the statistical distri- 
bution of the observable 8 is discussed. It can be written 



as 8 



2qxy + {q^ — l)y'^ = f{x,y) in terms of the 



matrix elements x, y and the mean value q of the Fano 
parameter. The probability distribution of 8 is given by 



P(8) = / dxP^^ix) / dyP,^{y)d{e-fix,y)) . (10) 



Introducing the scaled area 8 
^-1074+^2 and 6= [s + s 



tribution P(Q) can be written as 
1 



= a8 with a = 
"i(g2 - 1)] /4 the dis- 



P(8) 



exp(M8)i^o(|8|) 



(11) 



Additionally the 
has been intro- 



where Kq is the MacDonald function 
asymmetry parameter fj, = b[j + b^^ 
duced which determines how much P(8) deviates from 
a symmetric distribution. P(8) diverges logarithmically 
as 8 approaches zero. 

Fig. 1^ demonstrates the dependence of P(8) at fixed 
g = in the three characteristic regimes n < 0, fi — 
and /i > by varying s. For fi — the probability 
distribution is symmetric with respect to 8 = (solid 
line) . For /i < the area distribution is asymmetric with 
the flat side at negative values of 8 (dashed line). The 
opposite holds for > (dot-dashed line). 

As fi depends both on s and q a contour plot of fi as 
a function of {q,l/s) is presented in Fig. |[ The vari- 
able 1/s instead of s has been chosen for matters of dis- 
cussion as 1/s oo corresponds to the limit of strong 
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FIG. 4: Asymmetry parameter /i of the area distribution, 
Eq. (11) as a function of the average shape parameter q and 
s~^. The white dashed line marks the values (g, 1/s) where 
fj. = 0. 



motion in the closed channel is chaotic. This situation is 
frequently encountered in photodissociation of highly ex- 
cited rovibrational states of molecules or photoabsorption 
of atoms to auto-ionizing resonances. The closed channel 
subspace is modeled by random matrix theory. The dis- 
tribution of the Fano parameter q was derived for both 
the GOE and the GUE cases. The peak height distribu- 
tion and the resonance area distribution were discussed 
in detail for the GOE. 
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coupling between the closed and the continuum channel. 
The range of/iis— 1 < /i< 1. The white dashed line 
marks the subspace of parameters where ^^(6) is sym- 
metric (/i = 0), given by -I- = 1. If \q\ > 1 the 
asymmetry parameter /i is always positive regardless of 
1/s. For fixed \q\ < 1 and weak coupling (1/s small) the 
asymmetry parameter fj, is positive, too. If the coupling 
is enlarged fj, becomes negative above a certain value of 
1/s which depends on q. 

The expectation value of the scaled area distribution is 
given by (O) = /i (1 + 46^). Thus (9) is always positive 
when \q\ > 1. Note that the area under an individual 
resonance is positive if for its Fano parameter \q\ > 1 
holds. When \q\ < 1 the area associated with an indi- 
vidual resonance is negative. In contrast (0) can either 
be positive or negative for \q\ < 1 depending on the cou- 
pling strength to the continuum. Strong coupling to the 
continuum (1/s ^ oo) results in a negative value of (Q). 
In the limit fj, +1 Brcit-Wigner resonances dominate 
the spectrum. For fixed \q\ > 1 this limit is reached 
both for weak and large coupling to the continuum. For 
\q\ < 1 the Breit-Wigner limit is reached only for weak 
coupHng (1/s 0). For \q\ < 1 and strong coupHng 
to the continuum (1/s ^ oo) the expectation value of 
the area distribution is negative and window resonances 
with negative area dominate the spectrum in the limit 
/i^-1. 

Summary — Statistical properties of Fano resonances 
were derived for the situation of competing direct and 
indirect photoabsorption on the basis that the classical 
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